Influence of Cooper pairing on the inelastic processes in a gas of Fermi atoms. 
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Abstract 

Correlation properties in ultracold Fermi gas with negative scattering length and its impact on the three-body recombination 
. . . is analyzed. We find that Cooper pairing enhances the recombination rate in contrast to the decrease of this rate accompanying 
■ Bose-Einstein condensation in a Bose gas. This trend is characteristic for all interval of temperatures T<Tc. 
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The investigations of ultracold atomic gases give unique possibilities for studying quantum correlations and their 
role in the kinetics and dynamical properties of macroscopic systems. At early stage of these investigations it has 
been predicted that the rate of inelasticprocesses drastically decreases with the Bose condensation due to principal 
reconstruction of quantum correlations The local three-particle correlator K^, responsible for the recombination 
in a dilute gas, reduces with decreasing the temperature T < Tc and the ratio K3 (T) /A'3 (Tc) becomes equal to 
^1/6 as T — > 0. As the gas parameter na^ increases the effect decreases and practically vanishes at the liquid densities 
(a is the scattering length). 

For the first time, the effect has been observed in work Q. The authors studied the decay kinetics of a ultracold 
^'^Rb gas in a trap under the BEC realization. The experimental results turned out in a qualitative and quantitative 
agreement with the theoretical predictions. In fact, these experimental results were ones of first evidences for the 
quantum correlation formation in the process of the BEC kinetics with a limited lifetime of the system. 

The analogous phenomenon takes place in the 2D case at finite temperatures T < Tq when the condensate is absent 
This result reflects a key role of the local correlations which are close to the case of a genuine condensate. 

The condensate of Cooper pairs occurs in a two-component atomic Fermi gas with the attraction between particles 
at sufficiently low temperatures. It is interesting to reveal how the Cooper pairing and pair condensate affect the 
three-particle recombination rate. Analyzing this problem, we suppose that the attraction is a result of the Feshbah 
resonance for s-wave scattering particles. Herewith, the interaction is supposed to have nonzero value for the particles 
of different components only. 

The results obtained in the present work lead to a priory not evident conclusion, namely. Cooper pairing results in 
enhancing the probability of three-particle recombination. To make the investigation more transparent, we consider 
the low density gas, assuming the following conditions 



\ a \ kp « 1 and rokp « 1 (1) 
where tq is the characteristic spatial size of the interparticle interaction. 

In the case of negative scattering length a < the weakly bounded dimers, typical for a > 0, are not created. 
The recombination in this case is accompanied by the formation of a molecule at a deep energy level. Completely 
three-particle character of the relaxation in the case of the negative scattering length a < 0, including unitary region, 
is demonstrated in the work Q . The formation of the weakly bounded dimers in the case of the positive scattering 
length a > changes the kinetics of the strongly bound molecule formation essentially Q . The large released energy 
transfers to the kinetic energy of the molecule as a whole and to the third particle involved into the recombination 
process. Here the picture is analogous for Bose and Fermi particles and the transition rate is proportional to the 
three-particle correlator in the both cases. However, this correlator for Fermi particles demonstrates the increase with 
the Cooper pairing in contrast to its decrease in the BEC case. 

2. Considering three-particle recombination in a low density Fermi gas with the creation of a molecule in a strongly 
bound state, we present the Hamiltonian of the system in the form 

H^Ho + H' (2) 

where Hq is the Hamiltonian corresponding to elastic processes and H' is the Hamiltonian corresponding to inelastic 
recombination processes, respectively. 

In this expression ip^ is the creation operator of a molecule 




where ifm is the wave function of a molecule in its center-of-inertia frame. The Fermi operators of particle annihi- 
lation have the standard form 
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Conserving the pairwise structure of the interaction, we present V in the form V — U (l^i — 1^2) + U (1^3 
From (4) it follows that | T^i — 1^2 \ in (3) has a scale of the molecule size r*. The characteristic value of | 7^3 
in (3) is close to r* since in this case only the third particle should get the momentum compared in the magnitude 
with the momentum of the created molecule. Analyzing (3) in the frame which origin is at the center-of-inertia of the 
three-particle ensemble, we introduce the variables 



7^ = A 



^2, 



7^' = ^3-^2 



(6) 



The bare vertex for a product of three Fermi operators on the right hand-side with taking (6), (1) and the trans- 
position symmetry into account can be transformed as 
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(7) 



exp {i k i~p + i k 3~p'^ 1 — exp ^— i ^ fc 1 — fc 3^ (7^ — ~p 
^ exp (i"^!^ + it 3-^'^ (ti - 1^3) - -f') 



When the energy of a created molecule » ep, the momentum of the third particle and the momentum of the 
molecule are practically the same in the absolute magnitude. Taking into account the relation 
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for the vertex in (3), we have 
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Considering the Fermi gas with the large magnitude of negative scattering length | a |>> tq, and conserving 
limitation (1), we face with a sharp increase of vertex F under conditions of the Feshbah resonance. In fact, for 
the singlet pairs a quasi resonance state is realized in the continuous spectrum for energy e — > 0, whereas the real 
dimer bound state is absent. The solution of the Schroedinger equation for the two-particle problem demonstrates 
that, for I ~^ |< ro, the wave function obtains an additional large factor | a \ /rg. The probability amplitude to 
find simultaneously three fermions in the volume ^ requires a factor (| a | /rp) as compared with the case of 
noninteracting gas (c.p. [1], 0], Q)- Note that r» < tq. Correspondingly, vertex (8) obtains effectively an additional 
factor ~ (I a I /ro) . 

The integral in (8) can be represented as a sum of two integrals. Each of them is factorized as a product of two 
integrals over ~p and ~p' . Integrating in every case by parts and taking into account that | fc ^ |<< | |, we find. 



F = F(^; 



(9) 
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Here ^ is the numerical coefficient. In all cases we suppose the spherical symmetry of the functions U ("r^), <^m {'^)- 
As a result, F ("^) depends on the absolute value of vector 1^ alone. 

Considering the recombination transitions and, correspondingly, H' as a small perturbation, the number of transi- 
tions per unit time takes the form [h = 1) 

W = 2nJ2pi I H'fi f 5 {Ef - Ei) = j dt < H' (0) H' (t) > (10) 

H' {t) = exp (^iHot^ H' exp (^-iHotj 

Here the operator pi is the equilibrium density matrix defined by the Hamiltonian Hq. Summing over index f in 
(10), we use the inequality sp « Em- Then the block 

T -oo 



can be singled out in (10). Here the value | (g /m) is the overall kinetic energy of the created molecule and fast 
atom with the momenta equal in magnitude and opposite in the direction. The direct calculation of this expression 
gives 



^,,,,.((t-t.)(^'.-^0) 



As a result, for the probability (10) we find 



i=l,2,3 

where K^^^ is the three-particle correlator 

if (3) (ti, t^) =< (0) (0) (0)cfe3t (0) Cfe,; (0) Cfe,t (0) > (12) 
The value of the constant B is determined by the expression 

B=-^mq,\fiq,)\^ (13) 

The characteristic transition time 1/Em is small compared with the characteristic correlation time in (11), and it 
determines effectively the correlator (12) for the coincident time moments. 

3. We begin with the calculation of the probability W (11) at zero temperature T = 0. In the lack of interparticle 
interaction the correlator K^^^ reads 

\ ' */ it k2i featV k'i'^^i k3\k'^ feilfeg k3;k[J k^'tk'.^ 

Substituting this expression into (11), we find 

^o = B^|^n3, (14) 
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where n is the particle density with the fixed spin projection. If the interaction between particles is taken into 
account, correlator (12) should be averaged over the changed ground state. We start from the case when the Cooper 
pairing is absent and then find the correlator K'-'^^ in first order in the parameter | a | fci? << 1- 

The interaction Hamiltonian for low energy particles can be written as 

Hint = — V (l^)'cty ci Ct^ , C7+ 

— > — > 112* 
p 1, p 2 

where ~^ = — ~^ i and an additional condition for summing in this expression is ~^ i + = + "^2- In first 
order in the interaction the wave function of the ground state can be represented as 



The direct calculation of the correlator JiT^^-' with taking this reconstruction into account results in the appearance 
of four terms. If the factor fc 3 ) is involved, these terms give the identical summands in Eq. (11), 



and the correlator can be written as 



iff) {ta')--v{t 



K , k „ k I k 2 k -, fcg k I k 2 



n-^ .S-^ -r>, (15) 



where k = k 2 — k 2- this expression the term with a product of four functions n-^ in the numerator gives zero 
contribution to (11) due to the symmetry and can be omitted. 

We use the usual method of substituting genuine potential V (T") for an effective potential V ("?") with the same 
magnitude of the scattering length a for small energies e — > and conserve a possibility to employ the perturbation 
theory. The potential V {'f') can be taken as a simple spherical rectangular well with the depth Vq >> sp and the 
radius Rq. Applicability of the perturbation theory implies the validity of the inequality xRq << 1 where k = \/ mVo, 
herewith, the bound state in the well is absent. It can easily be revealed that, for these requirements, the relation 
Ro = 7] \ a \ with 77 >> 1 takes place as well as combined inequalities kp \ a \« kpRo « 1. As a consequence, the 

Fourier component obeys the equality V (^k — > 0^ « (47r/m) a = g and the function V (^k^ begins drastically to 
decrease for kRo ^1 {k >> kp)- The last point is important since with the substitution (15) into (11) and under 
the integration over momentum k the divergent of the second term in the brackets in (15) is eliminated when the 
dependence V (^k^ is taken into account. The magnitude of the integral proves to be ~ l/i?o and the ratio g/Ro 
becomes independent of the scattering length. In fact, the contribution of the corresponding term in the correlator 
K^^") determines some renormalization of the probability Wo (14). 

In the first term in (15) all the momenta obey inequality ki < kp and the effective interaction can be approximated 

as V (^k^ ~ Substituting iff-* into (11), we find the probability of the transition Wi and, at the same time, the 
general expression for W in the normal phase in the linear approximation in interaction 

Wn^Wn + Wi^Wo(^l-^ (11 - 21n2) kp \ a \^ (16) 

4. Consider now the contribution of the Cooper pairing to the formation of the three-particle correlator (12). For 
this purpose, averaging over the ground state, we should include the anomalous averages 

i4:f) = -<2i xd-r^ , > n-r^ s-r^ -^,s^, -^,6^ -f +■■■ (17) 

In this expression there are additional three terms resulting from the commutations of the momenta and giving the 
identical contribution to (11). 

Using the standard method of the u-v transformation for the calculation of the anomalous averages, see e.g. 0, we 
have at zero temperature 
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Here E}^^ = ^| A ^ + — and A is the Cooper gap in the one-particle excitation spectrum. Within the 
framework of the method considered the integral in the expression obtained is well known 



d^h A 2A 



(27r)' Ek, 



One can sec that the result holds for T ^ if the gap A is the temperature-dependent quantity A (T). As a result, 
after substituting Kg (17) into (11), we find directly 



From Eq. (18) it follows the main qualitative result: the Cooper pairing or condensate of Cooper pairs intensifies the 
three-particle recombination and thus the relaxation rate of the attractive Fermi gas. In this aspect the result differs in 
kind from the BEC case in a dilute atomic Bose gas in which, on the contrary, the rate of three-particle recombination 
decreases. Due to the behaviour of A as the function of the parameter kpa the three-particle recombination rate rises 
exponentially when kp \ a \ increases, at least up to kp \ a |< 1. The exterpolation of the results obtained from the 
region fc^ | a |<< 1 to the region kp | a |^ 1 demonstrates the real possibility of the experimental study of the effect. 

The leading temperature correction for Wn is connected with Wq and equal to (tt^/S) [T/ep) . Since — 
A(T = 0), it is readily seen that, for the interval T < T^, the temperature behavior of W is governed by the de- 
pendence Ws (T) (18). 
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